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TWO-DIMENSIONAL  UNSTEADY  LIFT  PROBLEMS 
IN  SUPERSONIC  FLIGHT 
By  ffex.  A.  Heaslet  and  Harvard  Lomax 


SUMMARY 

The  variation  of  preestire  distrihution  is  calculated  for  a  two- 
dimensional  supersonic  airfoil  either  experiencing  a  sudden  angle— of— 
attack  change  or  entering  a  sharp— edged  gust.  From  these  pressure 
distributions  the  indicial  lift  functions  applicable  to  unsteady 
lift  problems  are  determined  for  the  two  cases.  A  close  similarity 
is  shown  to  exist  between  the  calculated  functions  for  varying  free- 
stream  Mach  number  and  the  corresponding  functions  in  three- 
dimensional  incompressible  flow  for  varying  aspect  ratio. 

Results  eire  presented  which  permit  the  determination  of 
maximum  Increment  in  lift  coefficient  attained  by  an  unrestrained 
airfoil  during  its  flight  through  a  gust.  As  an  application  of 
these  results,  the  minimum  altitude  for  safe  flight  throu^  a 
specific  gust  is  calculated  for  a  particular  supersonic  wing  of 
given  strength  and  wing  loading. 


INTRODUCTION 

The  study  of  the  unsteady  lift  of  wings  in  an  incompressible 
medium  has  been  developed  along  two  different  lines.  In  refer¬ 
ence  1,  R.  T.  Jones  introduced  the  concept  of  indicial  lift 
functions  for  wings  of  finite  aspect  ratio  and,  using  as  a  basis 
the  work  of  Wagner  (reference  2)  on  the  two-dimensional  potential 
theory  of  airfoils  in  nonuniform  motion,  has  shown  how  the 
calculation  of  lift  under  various  conditions  of  motion  can  be 
effected.  In  reference  3,  Theodorsen  considered  nonsteady  motion 
in  iis  relation  to  the  general  theory  of  aerodynamic  instability 
and  the  determination  of  the  aerodynamic  forces  on  harmonically 
oscillating  airfoils.  This  latter  approach  has  been  extended  to 
Include  high-speed  problems  and  in  two  recent  papers  Garrick  and 
Rublnovf  (references  and  5)  have  given  results  on  flutter  and 
oscillating  air-force  calculations  for  wings  in  supersonic  flow. 
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The  presort  report  employs  the  method  of  attack  introduced  hy 
Jones  and  considers  the  case  of  a  two-dimensional  airfoil  moving 
supersonically  in  an  arbitrary  manner^  provided  the  assumptions  of 
small  perturbation  theory  are  satisfied.  The  principal  contribution 
lies  in  the  determination  of  Indicial  pressure  distributions  which 
are  readily  calculated  in  supersonic  motion  and  from  which  Indicial 
lift,  drag,^  and  pitching  moments  may  be  computed.  From  these 
results  the  Indicial  lift  functions  are  calculated  explicitly. 

The  methods  used  to  find  the  pressure  distributions  also  afford 
considerable  insight  into  the  same  problem  for  airfoils  at  subsonic 
speeds. 

As  an  application  of  the  analysis,  the  results  are  applied  to 
the  special  case  of  an  unrestrained  airfoil  entering  a  sharx>-edged 
gust.  The  resultant  forces  are  found  to  be  comparable  in  magni¬ 
tude,  for  Mach  numbers  in  the  neighborhood  of  1*3^  to  those  given 
in  reference  1  for  subsonic  incompressible  flow. 


s™bols 

a  speed  of  sound 

c  chord  length 

oi  section  lift  coefficient 

ci^  indicial  lift  coefficient  for  angle-of-attack  change 
Cl  indicial  lift  coefficient  for  wing  entering  gust 

g  acceleration  of  gravity 

M  free— stream  Mach  number 

m  mass 

Ap  difference  in  pressures  between  lower  and  upper  surfaces  of 
airfoil 

a  (fPoVo^) 

s  distance  measured  in  half— chord  lengths 
S  area  of  wing 

p  perturbation  static  pressure 
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time  In  seconda 

t  transformed  time  variable  (See  e(iixatlons  (7)») 

T  t/c 

u  perturbation  velocity  component  In  x  direction 

free-stream  velocity 

V  pertiirbatlon  velocity  component  In  z  direction 

Vq  z  component  of  velocity  of  gust 

W  weight  of  wing 

XyZ  Cartesian  coordinates 

a  angle  of  attack 

[X  2m/pMSc 

p  perturbation  density 

pQ  free— stream  density 

cr  area  over  which  surface  Integral  Is  evaluated 

$  perturbation  velocity  potential 


Subscripts 

o  free— stream  conditions 

1  variable  of  Integration 

u  upper  surface 

I  lower  surface 

AMLYSIS 

Derivation  of  Basic  Differential  Equation 

The  pressure  distribution  over  an  airfoil  In  a  compressible 
medium  3  s  obtainable  from  the  solution  of  a  boundary  value  problem 
associated  with  a  particular  second— order  partial  differential 
equation.  The  derivation  of  the  linearized  form  of  th3s  equation, 

obtained  under  the  assimiptlons  of  small  perturbation  theory, 
proceeds  as  follows : 
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Let  "be  perturbation  velocity  components  parallel, 

respectively,  to  the  Cartesian  axes  x*,  z*  and  denote  by 

p  perturbation  pressiire,  by  p  perturbation  density,  and  ^ 

by  a  the  velocity  of  sound.  Then  If  t*  denotes  time, 

Vq  13  the  constant  free -stream  velocity,  and  is  the 

constant  free— stream  density,  the  linearized  Eulei*lan  eqxiations 
are 


^  du  1  Sp 

St»  ^  Sx«  p  Sx* 

o 


(1) 


The  linearized  equations  of 


dv  1  Sp 

ox*  p^  oz* 

continuity  and  state  are,  respectively. 


St*  hx* 


(2) 


(3) 


After  the  introduction  of  the  perturbation  velocity  potential  $ 
integration  of  equations  (l)  yields 


5 $  ^  y  3$ 
at‘  °  bx' 


—  p  +  const. 

Po 


while  equations  (2)  and 


(3)  give  the  expression 


Sx'/  Poao^ 


(4) 


(5) 
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The  camhlnation  of  equations  (4)  and  (5)  leads  to  the  desired 
partial  differential  equation 


(1~M^)  ^  _ L 

ao  5z*dt*  ao^  dt*^ 


(6) 


where  M  is  the  free-^stream  Mach  numher. 

Equation  (6)  can  he  reduced  to  the  nonrialized  form  of  the  two- 
dimensional  wave  equation  of  mathematical  physics  hy  means  of  the 
transformation 


X  =  X*  —  Ma^t* 

z  =  z  *  > 

t  =  a  t’ 
o 


(7) 


In  these  variable s  the  equation  is  written 

^  ^  Q  /g\ 

In  accordance  with  the  assumptions  underlying  the  derivation  of 
equation  (8),  its  application  to  prohlems  in  airfoil  theory  is, 
of  course,  limited  to  cases  where  the  Induced  velocities  are  small 
compared  to  the  free— stream  velocity  and  the  effects  of  viscosity 
do  not  alter  the  results  of  the  potential  flow  solution. 

The  rectan^lar  coordinate  system  associated  with  equation  (6) 
is  fixed  in  the  wing  which  is,  in  turn.  Immersed  in  a  free  stream 
of  velocity  directed  along  the  positive  x— axis.  The  transforms^ 

tlon  introduced  in  equation  (7)  fixes  the  x,z  cordinate  system 
in  space  so  that  the  airfoil  moves  in  the  negative  x  direction  and 
the  free— stream  velocity  is  zero.  A  distortion  of  the  time  axes 
is  also  involved  so  that  the  differential  equation  appears  in 
canonical  form.  Equations  (6)  and  (8)  are,  of  course,  well  known 
in  the  theory  of  unsteady  motion  and  in  the  study  of  sound  waves. 

It  is  natural  that  these  two  fields  of  study  should  yield  the  same 
fundamental  equations  since,  in  small  perturbation  theory,  the  wing 
may  be  thought  of  as  a  distribution  of  acoustic  radiators. 
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Solution  For  Given  Boundary  Conditions 

The  bcundary  conditions  which  are  to  he  satisfied  have  the  saine 
property  as  those  encountered  In  steady— state  thlnr-alrfoll  theory j 
that  Is,  the  prescrlhed  data  are  given  In  the  z  =  0  plane.  The 
particular  prohlems  with  which  this  report  deals  are  those  of 
finding  pressure  dis  trihut  Iona  over  a  flat  plate.  Thus,  w  will 
he  specified  over  a  portion  of  the  z  =  0  plane  and,  elsewhere 
throughout  the  plane,  loading  must  he  zero. 

Two  houndary— value  problems  axe  to  he  considered;  first,  the 
case  of  an  airfoil  either  starting  from  rest  at  a  given  angle  of 
attack  a  or  experiencing  along  the  entire  chord  a  change  a  in 
stream  direction  without  a  pitching  motion,  and  second,  the  case 
of  a  constrained  wing  entering  a  sharp-edged  gust  with  a  vertical 
velocity  Wq.  In  the  former  case  the  motion  is  that  of  an  airfoil 
suddenly  sinking  without  rotation.  These  boundary  conditions  are 
more  readily  pictured  with  the  aid  of  figure  1.  Figure  1(a)  shows 
the  conditions  which  must  he  satisfied  in  order  to  solve  the  angle— 
of— attack  problem.  The  trace  of  the  leading  edge  of  the  wing 
traverses  the  line  x  =  -  Mt,  while  the  trailing  edge  lies  on 
X  =  c  —  Mt  where  c  is  chord  length.  The  region  hounded  by  these 
lines  and  the  line  t  =  0  is  the  region  occupied  by  the  airfoil 
as  time  passes.  Since  the  axes  are  fixed  and  the  airfoil  moves  in 
the  negative  x  direction,  the  velocity  at  which  the  airfoil 
travels  determines  the  inclination  of  the  loci  of  the  leading  and 
tralllng-edge  traces.  Over  the  "area”  occupied  by  the  airfoil  in 
the  x,t  plane,  w  must  equal  — Voa  and  elsewhere  no  jump  in 
pressure  can  occur.  The  gust  problem  (fig.  1(h) )  does  not  differ 
essentially  from  the  previous  problem  except  that  here  the  region 
over  which  the  modification  of  w  is  effective  is  not  entirely 
the  region  occupied  by  the  airfoil  hut  rather  the  region  occupied 
simultaneously  by  the  airfoil  and  the  gust.  Fixing,  for  convenience, 
the  edge  of  the  gust  along  the  h-axis,  this  axis  will  form  the  rlghtr- 
hand  boundary  of  the  region  over  which  w  =  —  Wq. 

Since  the  partial  differential  equation  is  linear  and  the 
solutions  are  therefore  additive,  these  boundary  conditions  clearly 
should  fit  the  following  physical  event:  a  wing  of  trapezoidal 
plan  form,  indicated  in  figure  2,  flies  at  a  steady  lift  and  angle 
of  attack  prior  to  t  =  0;  at  t  =  0  the  wing  either  experiences 
a  change  in  angle  of  attack  a  with  no  pitching  motion  or  enters 
a  sharp-edged  gust  of  constant  vertical  velocity  Wq,  the  gust 
extending  from  its  edge  to  all  negative  values  of  x.  The  wing  in 
each  case  is  then  restrained  so  that,  relative  to  the  original  wind 
vector,  the  wing  remains  at  an  angle  of  attack  a  or,  in  the  gust 
case,  continues  fixed  at  the  same  angle  of  attack. 
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The  solution  to  similar  houndary— value  prohlems  has  heen 
discussed  at  length  In  reference  6.  In  that  report  the  develop¬ 
ment  was  adapted  to  the  case  where  equation  (8)  represented  the 
steady— state  equation  for  wing  prohlems  in  three  dimensions  and 
the  chara-cteristic  cones  of  the  equation  had  the  immediate  physical 
interpretation  of  Mach  cones  or  infinitesimal  shock  disturbances. 
Greenes  theorem  was  applied  to  solve  the  boundary-value  problems 
involved  and  it  was  shown  that  the  solution  obtained  could  be 
Interpreted  as  a  surface  distribution  of  sources  and  doublets. 
Finally^  the  difficulties  arising  in  the  discussion  of  the 
singularities  on  the  characteristic  coiies  and  the  Integration  of 
the  supersonic  doublet  were  overcome  through  the  introduction  of 
an  Integration  technique  which  Involved  using  the  "finite  part" 
of  the  given  integrals. 

The  methods  and  conclusions  obtained  in  reference  6  can  be 
adapted  immediately  to  the  problems  discussed  herein  for  the 
mathematical  reasoning  remains  almost  identical.  The  physical 
interpretations  of  the  two  cases  must^  of  course,  be  modified. 

Thus,  the  characteristic  cones  of  the  differential  equations, 
traces  of  which  are  shown  in  figures  1  and  3,  ar©  no  longer  the  well- 
known  Mach  cones;  rather,  they  represent  the  distance  to  which  a 
disturbance  occurring  at  a  point  fixed  by  the  apex  of  the  cone 
will  travel  in  the  time  t.  Despite  the  fact  that  such  physical 
interpretations  are  undeniably  useful  in  understanding  and  applying 
the  results,  the  solution  of  the  basic  differential  equation  for 
the  boundary  values  Involved  is  quite  Independent  of  these  material 
dissimilarities.  It  follows  that  the  methods  developed  and  discussed 
in  detail  in  reference  6  can  be  applied  directly  to  the  given 
problems  with  only  minor  changes  in  notation. 

Eeferring  to  figure  1,  it  is  evident  that  for  supersonic  flight 
the  air  ahead  of  the  wing  is  unaffected  by  the  approach  of  the  wing 
and,  further,  that  the  induced  velocities  on  the  upper  surface  of 
the  airfoil  are  independent  of  the  shape  of  the  lower  surface. 
Consequently  the  pressure  distribution  will  be  found  on  the  upper 
surface,  as  if  the  airfoil  section  were  symmetrical,  and  then,  for 
the  flat  plate,  the  pressure  distribution  on  the  lower  surface  will 
be  equal  in  magnitude  and  opposite  in  sign.  As  in  reference  6,  the 
solution  to  such  a  problem  can  be  obtained  from  a  distribution  of 
sources. 

In  the  actual  computation  of  the  pressures  over  a  section 
traveling  at  supersonic  speeds,  certain  regions  are  conveniently 
defined.  These  regions  depend  on  the  relative  slope  of  the  traces 
of  the  leading  and  trailing  edges  and  the  trace  of  the  character¬ 
istic  cone  in  the  x,t  plane  (fig.  3).  The  perturbation  velocity 
potential  is  given  by  the  formula 
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dxi  dti 

«/  (t-ti)2  -  (x-Ii)2  -  z2 


(9) 


•where  a  is  'the  area  in  the  x,t  plane  of  the  region  occupied  by 
■the  wing  section  and  bounded  by  the  trace  of  the  forecone  from  the 
point  i,z,t  ..  If  p^  and  p^  denote,  respectively,  pressures  on 

the  upper  and  lower  surface  of  the  wing 


Ap 

4 


p  — p  ii-  8$  \ 

iLiu  =  — (  —  +  v„  —  ) 


4  8$ 

VqM  8t 


(10) 


Using  eguations  (9)  and  (10),  direct  calculation  shows  that  for  a 
sudden  angle— of-attack  change,  the  following  relations  hold 

Region  A  (be'tween  lines  2  =  —  Mt,  x  =  —  t,  and  x  =  c  —  Mt) 


Ap  _  ka 
^  1 


(11a) 


It  follows  that  in  this  region  the  steady-state  Ackeret-type  load 
distribution  has  been  at'tained. 

Region  B  (between  lines  x  =  —  t,  x=t,  and  x  =  c  —  Mt) 


Ap  Ikt  r  1 

^  yM2-l 


Mx+t 

arc  cos - + 

x+Mt 


kM 


+  arc  sin  J  (11^) 


Region  C  (be'tween  lines  x  =  t. 


t  =  0,  and  X  =  0  —  Mt) 


Ap  _  4a 
g  M 


(11c) 


The  result  ob-tained  for  Region  C  is  of  particular  in-terest  since 
it  holds  for  airfoils  at  subsonic  as  well  as  supersonic  speeds. 
Moreover,  the  mechanics  of  the  interaction  between  the  airfoil  and 
the  fluid  are  such  that  other  methods  of  derivation,  furnishing 

added  insight  into  the  nature  of  the  phenomenon,  may  be  developed. 

Consider  a  flat  plate  of  infinite  aspect  ratio  flying  at  a  velocity 
Vq  either  greater  than  or  less  than  the  velocity  of  sound  a^  in 
the  undisturbed  air.  The  airfoil  is  assumed  to  undergo  a  change  in 
its  motion  at  the  time  t  =  0  so  that  subseguent  to  this  time  it 

has  increased  its  angle  of  attack  by  the  amount  a.  It  follows  that 

the  sudden  increment  in  lift  can  be  calculated  from  a  knowledge  of 
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the  Irduced  effects  on  the  air  produced  hy  an  added  vertical  velocity 
of  the  plate  eq,ual  to  — Voa.  As  a  result  of  this  vertical  Kotlon  two 
plane  Eaylelgh  waves  will  emanate  from  the  plate,  a  compresslcn  wave 
from  the  lower  surface  and  an  expansion  wave  from  the  upper  surface. 
The  velocities  of  the  wave  fronts  are  equal  to  a^  while  the 
induced  velocities  in  the  waves  are  equal  in  magnitude  to  Vga. 

The  lift  on  the  wing  can  he  determined  in  two  ways:  from 
impulse  relations  and  from  energy  considerations.  In  the  former 
case,  assume  that  the  forces  per  unit  span  on  the  upper  and  lower 
surfaces  are  fu  and  fj,  respectively.  After  an  elapse  of  time 
At  the  wave  fronts  have  advanced  a  distance  aoAt  and  each 
includes  a  mass  of  air  equal  to  p^caQAt  per  unit  span.  From 
Newton's  second  law  of  motion 


(f ^-fy)At  =  (PoCa^At)  (2V'Qa) 


and,  converting  to  lift  coefficient, 

f^— fu  ^ 

Since  the  force  is  distributed  uniformly  along  the  chord  this  resuilt 
is  the  equivalent  of  equation  (11c) . 

In  the  development  of  the  theory  of  plane  waves  of  small 
amplitude  (see,  e.g.,  reference  7)  Eaylelgh  and  Iamb  have  shown 
that  the  energy  in  a  wave  is  divided  equally  into  kinetic  and 
potential  energy.  Denoting  kinetic  energy  by  T, 


w^dxdydt 


where  w  is  the  perturbation  velocity  within  the  wave.  Since  the 
energy  induced  in  the  wave  must  result  from  work  done  on  the  plate, 
it  can  be  seen  that 


=  I  Po  7  Vo^a^cdx 
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or 

(f2~f^)At  =  2PoVoa  ca^At 

This  equation  is  in  agreement  with  the  one  obtained  previously. 

As  a  consequence  of  equation  (11c)  it  follows  that  the 
starting  lift  coefficient  of  an  airfoil  Is  equal  to  ka/M  for  both 
subsonic  and  supersonic  flight.  The  magnitude  of  c^  thus 
increases  as  M  becomes  smaller  and  for  incompressible  theory^ 
where  the  velocity  of  soiind  is  indefinitely  large,  must  necessarily 
become  infinite.  This  fact  was  known  previously  along  with  the 
understanding  that  the  indicial  lift  function  experiences  an 
infinite  discontinuity  at  t  =  0.  (See  fig.  6.)  For  values  of  M 
other  than  zero  the  starting  lift  is  finite  and  a  continuous  lift 
function  results. 

For  a  sudden  gust  with  vertical  velocity  Wq  the  following 
e2i)ressione  can  be  found  for  the  corresponding  regions 

Eegion  A 


—  =  ^'*^0 

^  VqV'm^-I 


(12a) 


Region  B 

Ap  4wq  Mx+t 

—  =  - arc  cos  - 

Region  C 

Ap 

T  ” 


APPLICATIONS  AND  DISCUSSIONS 


Discussion  of  Load  Distributions 

Figure  k  shows  the  variation  of  the  loading  on  a  section  which, 
while  traveling  at  supersonic  speed,  is  suddenly  deflected  to  a  new 
angle  of  attack.  The  loading  varies  according  to  equation  (11)  for 
each  of  the  three  regions  A,  B,  and  C  of  figure  3*  At  t  =  0  the 
pressure  is  discontinuous.  Jumping  from  its  original  value.  Just 
before  the  sudden  deflection  in  angle  of  attack,  to  ka/u  Just 
after  the  deflection.  Figure  3  shows,  however,  that  the  initial 
l(Dad  distribution  is  modified  over  the  forward  portion  of  the  section 
as  time  increases  since  regions  A  and  B  must  be  considered. 


(I2t) 


(12c) 
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For  t  'the  ving  lies  entirely  in  region  A  so  that 

Ap  _  4a 
1  -/M^-l 

and  the  loading  has  attained  a  static  value  agreeing  with  the  Ackeret 
type  of  distribution .  The  loading  in  region  B  varies  between  the  two 
constant  values  of  regions  A  and  C,  dipping  below  that  of  region  C 
and  having  its  mininrum  value  at  x  =  0. 

The  loading  produced  upon  entering  a  sharp-edged  gust  is  pictured 
in  figure  5.  The  loading  in  region  C  is  zero, since  that  portion  of 
the  wing  is  unaware  of  the  change  in  stream  conditions.  Over  the 
forward  portion  of  the  airfoil  (region  A)  the  Ackeret  type  loading 
corresponding  to  the  modified  angle  of  attack  is  in  evidence  and  for 

t  >-^  extends  over  the  entire  chord  of  the  wing.  In  region  B 
M— 1 

th©  loading  experiences  a  reduction  in  magnitude  from  the  value  over 
the  forward  portion  of  the  wing* 

The  load  dlstrlhutlons  which  have  "been  developed  were  obtained 
for  flight  velocities  in  the  supersonic  regime.  It  is  apparent, 
however,  that  the  basic  differential  eq,uation  is  not  restricted  to 
the  case  where  M  >1  and  that  the  method  of  analysis  affords  a 
means  whereby  transient  load  effects  can  be  studied  for  subsonic 
speeds.  The  essential  difference  between  the  latter  problem  and 
the  results  derived  here  lies  in  the  relative  position  of  leading— 
and  trailing— edge  traces  in  the  x,t  plan©  and  the  trace  of  the 
characteristic  cone.  Thus,  for  subsonic  flight,  the  trace  x  =  — t 
does  not  cut  across  the  region  occupied  by  the  airfoil;  whereas  the 
cone  stemming  from  the  tralling-edge  point  c,0  does.  A  qualitative 
picture  of  the  problem  is  obtained  if  the  analogy  between  the  non- 
steady  two-dimensional  case  and  three-dimensional  wing  theory  is 
used.  The  loading  functions  given  in  equation  (ll)  are  equivalent 
to  loading  existing  on  a  swept-forward  tip  of  a  three-dimensional 
wing.  Thus,  in  figure  1(a),  x  can  represent  distance  measured 
spanwise,  t  can  represent  distance  measured  chordwise,  and  the 
shaded  area  can  represent  a  portion  of  the  plan  form  of  the  wing. 

Using  this  analogy,  the  loading  which  has  been  determined  is  merely 
load  distribution  over  the  swept-forward  tip  of  a  wing  with  constant 
chord  and  supersonic  leading  edge.  When  the  case  of  the  airfoil 
section  traveling  at  subsonic  speeds  is  to  be  considered,  the  problem 
becomes  one  of  determining  the  loading  over  the  swept-forward  tip  of 
a  wing  with  constant  chord  and  subsonic  leading  edge. 
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Development  of  Indiclal  Lift  Functions 
Since  section  lift  coefficient,  cj^,  is  given  ty  the  expression 

1  r°Ap 

°L  =  -  /  —  ^ 

0  Jo  q, 

the  relations  presented  in  equations  (11)  and  (12)  are  sufficient 
for  the  determination  of  cj^,  indicial  lift  coefficient  for  change 

in  angle  of  attack,  and  ct,  indiclal  lift  coefficient  for  an 

airfoil  entering  a  gust.  As  a  result  of  direct  integration,  the 
following  results  are  obtained. 


First  time  interval 


jia 

M 


CV, 


Second  time  interval  ° 


(13a) 

(Ika) 


°  <t<  ° 


~  r  —f'— 

Jt  L  MV2 


+  arc  sin 


1+M 


C““Mt 


M-1 


Vir-i 


arc  cos 


t+Mc-tM^ 


t2-(c-tM)2  j 


(131>) 


‘8 


=  f'l 

jtVoC  \2 


+  arc  sin 


!Ehird  time  interval 
ka 


c-Mt\  ^  ItWo  Mc+t®-M®t 


arc  cos 


(l4b) 


M-1 


(13c) 
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Values  of  the  lift  functions  are  plotted  In  figure  6  as 
functions  of  s,  the  distance  traveled  by  the  airfoil  measured  In 

half-chord  lengths  where  s  =  — .  The  curves  shown  were  calculated 

for  values  of  M  eq^ual  to  1.2,  1.31,  and  1.46,  since  the  asymptotic 
values  of  ci^  and  cl^  for  the  three  cases  agree  with  the  values 

given  In  reference  1  for  the  subsonic  wing  at  aspect  ratios  of  00, 

6,  and  3,  respectively.  No  direct  analogy,  of  course,  can  be  made 
between  the  two  cases.  It  Is,  however,  worthy  of  note  that  the 
variations  In  the  Indlclal  fimctlons  for  the  supersonic  case  are  of 
the  same  order  of  magnitude  as  those  found  In  the  finite-span 
Inc onLpre s s ilDle  case. 

Fram  a  Imovrledge  of  the  lift  function  resulting  from  a  sudden 
unit  angle  of  attaclr,  it  is  possible  to  express  the  lift  corresponding 
to  a  given  variable  motion  by  considering  the  given  motion  as  being 
composed  of  infinitesimal  steps  and  summing  the  lifts  corresponding 
to  each  step.  Mathematically,  the  problem  corresponds  to  the  use 
of  the  indlcial  admittance  in  determining  the  current  response  for 
an  electric  network  and  leads  to  the  so-called  superposition  theorem 
which  can  be  written  in  the  form 


cL(t») 


a-(o)  +y  CI^ 


(f-T«) 


da(  t'  ) 
dT» 


dr* 


(15) 


The  primes  on  the  variables  In  this  equation  Indicate  that  true  time 
is  used. 


Motion  of  Airfoil  in  Gust 


The  results  which  have  been  obtained  will  now  be  applied  to 
determine  the  forces  on  an  unrestrained  airfoil  entering  a  gust. 
Since  the  motion  of  the  wing  is  not  prescribed  it  becomes  necessary 
to  equate  the  dynamical  forces  in  order  to  relate  the  variables 
involved.  Neglecting  pitching  moment  and  using  Newton’s  second  law 
of  motion. 


m 


dw 

dF 


force 8 


(16) 


where  w  is  the  vertical  velocity  of  the  wing,  m  is  the  mass  of 
the  wing,  and  the  forces  to  be  summed  result  from  the  lift  on  the 
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wing  and  the  impressed  force  resulting  from  the  action  of  the  gust. 

By  means  of  equation  (I5),  eamtion  (I6)  can  be  rewritten  in  the  form 

(f-T*)a(T-)dT-  =  cig  (f)(iS  (17) 


Introducing  a  change  of  -variables  such  that 


and  setting 


_  t' 

T  =  —  ao. 


■w  _  „ 

—  =  a  , 

T 


_  T* 

Ti  =  —  ao 


H  = 


2m 


PqMSc 


the  equation  becomes,  finally. 


Vo  "g 


(18) 


Since  a(o)  =0,  ‘equation  (18)  can  be  integrated  to  give: 
•p  T 


a 


1  / 
Vo  d  Jo 


T  n  P  1 

CLg  (Ti)dTi  ^  j/  °La  ^  ° 

which  is  an  integral  equation  of  the  second  kind  with  a  variable 
upper  limit.  The  solution  to  this  equation  can  he  obtained  quite 
satisfactorily  by  means  of  Liouville’s  method  of  successive 
substitution.  Using  the  relation  c^  =  and  performing  the 

proper  manipulations  gives : 

T 

^  CL  -  («  -  U 

+  i  r^Cr  (T-Ti)dTir^CL  (T1-T2)  Cx  (T2)dT2 -  (20) 

d^^  Hi,  ® 

m  ^d^  M^— 1 

Equation  (20)  is  known  to  converge  uniformly  for  T< 

and  in  the  applications  of  this  report  the  maximum  lift  was  always 
experienced  in  the  region  of  convergence. 

The  values  of  lift  coefficient  Cl  determined  from  the  solution 
of  equation  (20)  are  shown  in  figures  7(a),  7(t),  and  7(c)  for 
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various  values  of  [x  and  for  M  =  1.2,  1.31j  and  1.46.  Figure  8 
shows  the  variation  of  the  mxlinum  lift  coefficient  attained  plotted 
as  a  function  of  the  density  parameter  \x  for  the  same  three  values 
of  M.  Figure  9  furnishes  a  coii5)arlson  "between  the  maximum  lift— 
coefficient  Increment  given  In  reference  1  for  an  aspect  ratio  of  6 
and  the  corresponding  value  calculated  In  the  present  report  for 

M  =  1.31*  The  results  are  plotted  as  functions  of  |ij  =  to 

correspond  with  the  density  parameter  used  "by  Jones.  The 
correspondence  which  was  noted  for  the  Indlclal  lift  functions  with 
Mach  number  replacing  wing  aspect  ratio  Is  still  In  evidence. 


Forces  Developed  on  G-lven  Wing 

As  an  example  of  the  uses  to  which  the  results  Just  obtained 
can  be  applied,  consider  an  airplane  with  wing  of  plan  form  such 
as  the  one  shown  in  figure  2.  Assume  a  wing  loading  of  40  poimds 
per  square  foot,  a  chord  length  of  8  feet,  and  let  It  be  specified 
that  the  wing  is  flying  at  a  Mach  number  equal  to  1.2  and  that  the 
wing  Is  built  to  withstand  forces  producing  accelerations  between 
—3  ajnid  5  times  gravitational  acceleration.  It  Is  proposed  to  find 
at  what  altitudes  the  wing  may  be  subjected  safely  to  a  gust 
possessing  a  vertical  velocity  of  50  feet  per  second. 

If  F  denotes  total  force  on  the  wing, 

F  =  W  +  Lg  (21) 

where  W  is  wing  weight  and  Lg  is  the  total  lift  produced  by  the 

gust.  If  ^  is  the  maximum  increment  in  lift  coefficient 

max 

attained  in  a  unit  gust  and  A  is  the  acceleration  factor  measured 
in  multiples  of  g,  then  equation  (21)  is  expressible  in  the  form 


A  =  1  + 


(-^l) 


max 


Wq  S 

—  —  0  V 
2  W  ^0  0 


Since  It  is  required  that 


A— l|  ^  h,  it  follows  that 


(^l) 


max 


6.k 

PoVo 


(22) 


Figure  10  shows  the  limit  curve  of  plotted  as  a  function 

of  flight  altitude.  From  a  knowledge  of  p  =  2W/pQgMSc,  however. 
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the  value  of  (^L)tnnY  actually  attained  by  the  wing  entering  the 

UIcxJL 

gust  can  te  calculated.  Such  values  are  also  Included  in  figure  10 
and  indicate  that^  under  the  given  conditions  the  wing  should  not 
fly  at  an  altitude  less  than  aj)prox3 irately  28,000  feet. 


Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Moffett  Field,  Calif. 
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(a)  Supersonic  wing  after  a  sudden  angle  of  attack 
change. 


X 


(b)  Supersonic  wing  entering  a  sharp  edged  gust 
at  x~  o 

Figure  /.  ~  Boundary  conditions  for  supersonic  wings. 
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I  H 

'^^Trace  of  Mach  cones 
\  from  leading  edge  tip 

\ 

\ 

\ 


Figure  2.-  Type  of  plan  form  studied  in  analysis. 


Figure  J  -  Sketch  indicating  locations  of  regions  A,  B  and 
C  and  time  int  ervals  used  in  analysis. 
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AP 

q 

Regions 

4  oc 

A 

B 

C 

yjM^-l 

4oc 

1 

M 

leading  edge  trailing  edge 

(a)  Time  interval  I,  0  <t/c  <  l/(M^I ) 


(b)  Time  interval  II,  l/(M^I)<  t/c  <I/(M-I) 


(c)  Time  interval  HI,  l/(M-l)<  f/c 

Figure  4.-Variations  of  Ap/q  with  chordwise  station  after  sud¬ 
den  change  in  angle  of  attack.  Regions  defined  in  figure  3. 
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(a)  Time  interval  I,  0  <t/c  <I/(M-^  !) 


(b)  Time  interval  l/(M+l)<  t/c  <  l/(M-l) 


M 

Q 


4w» 

^oyfM^I 


Region 

A 


leading  edge 


trailing  edge 


(c)  Time  interval  HI,  l/(M-  IJ  <  t/c 

Figure  5-Variation  of  Apy/q  with  chordwise  station  after  enter¬ 
ing  gust.  Regions  defined  in  figure  3. 
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Increment  in  lift  coefficient 
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Half  -  chords  traveled,  s 

Figure  7  -  Variatiation  of  increment  of  lift  coefficient  during 
flight  through  unit  sharp-edged  gust.  (a)M=l.20. 
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Figure  7.  -  Continued,  (b)  M- 1.31 


Increment  in  lift  coefficient,^  zlc. 
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Half -chords  traveted^s 


Figure  7.  -  Concluded,  (c)  M~l.46 
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Figure  9 .-Comparison  of  maximum  increments  in  lift  coef¬ 
ficients  during  flight  through  sharp-edged  gust  between 
subsonic  wing  and  supersonic  wing. 


Figure  10.-  Safe  altitudes  for  wing  with  loading  of  40  lb/ft‘ 
flying  at  M=I20,  based  on  50  ft/sec.  gust  velocity 
and  maximum  load  factors  of  5 g  and  -3  g. 


